Investigation of pulsed laser induced dewetting in nanoscopic metal films 
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Hydrodynamic pattern formation (PF) and dewetting resulting from pulsed laser induced melting 
of nanoscopic metal films have been used to create spatially ordered metal nanoparticle arrays with 
monomodal size distribution on SiCVSi substrates. PF was investigated for film thickness h < 
7 nm < laser absorption depth ~ 11 nm and different sets of laser parameters, including energy 
density E and the irradiation time, as measured by the number of pulses n. PF was only observed 
to occur for E > E m , where E m denotes the /i-dependent threshold energy required to melt the 
film. Even at such small length scales, theoretical predictions for E m obtained from a continuum- 
level lumped parameter heat transfer model for the film temperature, coupled with the 1-D transient 
heat equation for the substrate phase, were consistent with experimental observations provided that 
the thickness dependence of the reflectivity of the metal-substrate bilayer was incorporated into the 
analysis. The model also predicted that perturbations in h would result in intrinsic thermal gradients 
dT/dh whose magnitude and sign depend on h, with dT/dh > for h < h c and dT/dh < for 
h > h c ~ 9 nm. For the thickness range investigated here, the resulting thermocapillary effect was 
minimal since the thermal diffusion time th < the pulse time. Consequently, the spacing between 
the nanoparticles and the particle diameter were found to increase as h 2 and ft 5 / 3 respectively, which 
is consistent with the predictions of the thin film hydrodynamic (TFH) dewetting theory. PF was 
characterized by the appearance of discrete holes followed by bicontinuous or cellular patterns which 
finally consolidated into nanoparticles via capillary flow rather than via Rayleigh-like instabilities 
reported for low temperature dewetting of viscous liquids. This difference is attributed to the high 
capillary velocities of the liquid metal arising from its relatively large interfacial tension and low 
viscosity as well as the smaller length scales of the liquid bridges in the experiments. The predicted 
liquid phase lifetime tl was between 2 — 15 ns, which is much smaller than the dewetting time 
Try > 25 ns as predicted by the linear TFH theory. Therefore, dewetting required the application 
of multiple pulses. During the early stages of dewetting, the ripening rate, as measured by the rate 
of change of characteristic ordering length with respect to n, increased linearly with E due to the 
linear increase in t£ with increasing E as predicted by the thermal model. The final nanoparticle 
spacing was robust, independent of E and n, and only dependent on h due to the relatively weak 
temperature dependence of the thermophysical properties of the metal (Co). These results suggest 
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that fast thermal processing combined with the unique thermophysical parameters of metals can lead 
to novel pattern formation, including quenching of a wide range of length scales and morphologies. 

I. INTRODUCTION 



Physical phenomena that impose patterns with characteristic length scales can be utilized for the self- 
assembly of spatially ordered metal nanoparticles which have various applications, including in magnetics, 



catalysis, nanophotonics and plasmonics lllliyL 



fttDLIaLZD. Nonlinear dynamical instabilities which evolve 



into a robust stable state characterized by a narrow range of length scales are especially well-suited for 
this purpose. Typically such self-assembly of ordered patterns on surfaces and thin films arises from the 
competition between various physico-chemical effects. For example, surface rippling under ion irradiation 
can be attributed to an instability resulting from the competition between irradiation induced ion erosion 
and smoothing due to surface diffusion. The length and time scales of the ensuing patterns are determined 
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Alternately, 



by ion flux, surface temperature and surface diffusion parameters |& |9L 
nanopatterns could result from instabilities of a thin fluid layer leading to spatio-temporal re-organization 
of the fluid with well defined length and time scales that depend on the thermophysical material properties 
such as interfacial tension, fluid/substrate contact angle, fluid viscosity and, for ultrathin films, long range 
dispersion forces such as van der Waal's interactions. An example of pattern formation by such instabilities 
is dewetting of thin films caused by a hydrodynamic instability which occurs when attractive intermolecular 



forces overcome the stabilizing effect of interfacial tension B15L 



30, fl. 



Under such conditions, film 



thickness fluctuations are spontaneously amplified eventually leading to film rupture and the formation of 
drops/particles with well defined spatial order Jl^ll . 

Dewetting dynamics leading to particle formation have been studied in detail in polymer films that are 



in the liquid state close to room temperature 



BE, 



22, 



23Q . In comparison, studying dewetting in metals 



is challenging due to the high melting point of most metals. Metal dewetting under conventional long-time 
annealing is complicated because various surface processes including metal-substrate chemical interactions 



and metal diffusion into the substrate can influence dewetting 



techniques such as ion irradiation j26l 



25|]. On the other hand, rapid melting 



2711 or ns pulse laser heating could minimize such undesirable ef- 



fects and permit experiments over practical time scales |28fl. Approximately 10 years ago it was shown 



that single ns laser pulse meltingof various metals (Au, Cu, Ni) on a Cr metal layer resulted in dewet- 
ting patterns with spatial order |22, 30]. However, pattern formation and self-assembly of spatially ordered 
patterns under pulsed laser dewetting of metal films deposited directly on Si02 (a substrate of choice for 
many opto-electronic applications) have not been pursued in great detail, especially in nanoscopic metal 
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films with thickness h < the laser absorption depth a," 1 . Furthermore, pattern evolution under the ap- 
plication of thermally-independent multiple pulses (i.e. the film undergoes melting during the pulse and 
resolidification between pulses) with varying laser parameters to understand the dynamics of metal dewet- 
ting as well as the length scales of the variety of different patterns have not been investigated. This is of 
fundamental importance since the unique thermophysical properties of metallic liquids, including their high 
surface tension 7 of O(l) J/m 2 and large thermal conductivities, distinguish them from the typical liquid 
polymers (7 ~ 0(1CT 2 ) J/m 2 ) investigated in the past. Therefore, liquid metals could present novel pat- 
terning mechanisms and/or characteristics. Another important question of broader fundamental interest is 
related to the applicability of continuum-level conservation equations (for heat, mass and momentum) at the 
nanoscopic length scales (h < a" 1 ~ 11 nm) used in the present experiments. Note that the application 
of atomistic or mesoscopic techniques to problems which involve phase change and surface deformation 
under fast thermal cycles is not entirely straightforward and first principles-based atomistic/mesoscopic 
simulations of such phenomena are beyond present capabilities. Hence, the importance of developing ex- 
perimentally validated continuum-level models for such processes cannot be over-emphasized. 



Recently we showed that mu 



long-range spatial order [31 



32, 



tiple ns pulsed laser melting of nanoscopic Co films leads to short- and/or 



3311 . We provided preliminary evidence that dewetting by multiple ns pulse 
laser melting shows characteristics of a thin film hydrodynamic (TFH) instability, such as the spinodal 
process 03 lL This was based on the observation that the final nanoparticle morphology, quantitatively 
described by the nearest-neighbour (NN) particle spacing (R) and diameter (D), is consistent with the 
prediction of the linear TFH theory which was modified to incorporate long range inter-molecular attractive 
forces of van der Waal's type. However, the quantitative role of various laser parameters including the 
energy density E and time (as measured implicitly by the number of pulses n) on pattern evolution was not 
investigated in detail. In this work we present a detailed account of pattern formation and dewetting under 
multiple ns pulsed laser irradiation of Co nanofilms (h < 7nm < a" 1 ~ 11 nm) on Si02 in the light of 
first-principles thermal modeling and TFH theory (Sec. IIITl) and discuss novel physical mechanisms specific 
to dewetting of nanolayers of liquid metals. 

The remainder of the paper is organized as follows. In Sec. [IT]the experiment and modeling details are 
presented. In Sec. MI Al we show that by accounting appropriately for the changes in the reflectivity of 
the substrate-metal bilayer with respect to h, quantitative agreement can be obtained between experimen- 
tal observations and the predictions of continuum-level theory for the /j-dependence of the melting energy 



threshold, pattern length scale and particle size. This observation corroborates atomistic simulations of 



leat 



transport for simpler fluids under shear flow which, for liquid layers with several (6-7) atomic layers [34], 



yielded thermal profiles consistent with predictions of continuum-level analysis. For the metal films used in 
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this work, this translates to h < 1 nm. The thermal model also predicted the development of intrinsic thermal 
gradients dT/dh whose sign is positive for h < a critical value h c ~ 9 — 13 nm and negative for h > h c . 
This would imply that fhermocapiilary effects (i.e. flow induced by surface tension gradient along the per- 
turbed interface which samples an /z-dependent temperature field) could also be an important mechanism 
of pattern formation. For the experimental thickness regime investigated, the thermal diffusion time scales 
t~h < pulse time t p = 9 res and therefore the influence of dT/dh was observed to be minimal. However, 
since th oc h 4 , thermocapillary mechanisms could be important for thicker films, a hypothesis that needs 
to be verified by future experimental studies. In Sec. HUB I we discuss the quantitative behavior of pattern 
evolution as a function of E and re. Extremely fast thermal cycles (heating/cooling rates of O(10 10 K/s)) 
were predicted by the thermal modeling, implying that each laser pulse was thermally independent of one 
another (repetition rate: 50 Hz). Furthermore, the model predicted liquid lifetimes tl of 2 — 15 ns, while the 
linear TFH theory [15] predicted a dewetting time scale td > 25 ns implying that multiple pulses would be 
required to observe the instability, as confirmed by experiments. We suggest that once dewetting was ini- 
tiated by melting of the thin film, subsequent pulses continued to foster the instability and each laser pulse 
captured snapshots of its progression. Experimentally, the dewetting patterns were characterized by short 
range order (SRO) at all observed stages of dewetting with the patterns typically evolving with increas- 
ing number of laser pulses re from discrete holes to bicontinuous or cellular structures and eventually into 
nanoparticles. These nanoparticles form by capillary flow rather than a Rayleigh-like breakup typically ob- 
served in polymer liquids. This behavior can be attributed to the large capillary velocity (ratio of interfacial 
tension 7 to fluid viscosity n) for the metal nanolayers resulting from its unique thermophysical parameters, 
including a large 7 = 1.88 J/m 2 , small n = 4.45 x W~ 3 Pa — s and small pattern length scales due to 
its large Hamaker constant A ~ 1CT 18 J. This final nanoparticle state was robust, characterized by SRO 
in NN spacing R and a monomodal size distribution with diameter D with the observed trends of R ~ h 2 
and D ~ /1 5 / 3 respectively. This robustness could be explained by two effects: (i) the relatively weak tem- 
perature dependence of the thermophysical parameters of Co; and (ii) the minimal particle coarsening due 
to surface diffusion since the diffusion time scale is much greater than that of the thermal processing time 
t p . In Sec. MI CI we show that the pattern ripening rate L/n, where L is the characteristic pattern length, 
increased linearly with E at the early stages of patterning which could be explained by the predicted linear 
increase in tl with E. Overall, these results show that pattern formation and dewetting in liquid metal 
nanolayers under fast thermal cycles can result in novel mechanisms and, experimental observations in 
conjunction with theoretical analyses of thermal cycles and TFH instabilities offer quantitative guidelines 
for achieving desired quenched-in morphologies with predictable length scales. 
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II. PULSED LASER PROCESSING: EXPERIMENT AND THERMAL MODELLING 

A. Experimental details 

Cobalt metal films with thickness ranging from 1 < h < 7 nm were deposited at the rate of 1 nm / min 
onto optical quality SiC^/Si substrates at room temperature. The substrates were commercially obtained 
with a 400 nm thick thermally grown oxide layer on polished Si(100) wafers. Deposition was done under 
vacuum conditions (2 x 10~ 8 Torr) by e-beam evaporation from a TECTRA e" -flux mini electron-beam 
evaporator. The wafers were degreased in acetone, methanol, and DI water prior to film deposition. The film 
thickness was controlled through the deposition rate as measured using an in-situ Inficon XTM/2 deposition 
monitor. The final thickness was established by calibrating concentration measurements of Co film using 
energy dispersive X-ray spectrometry (EDS) in a scanning electron microscope (SEM) to measurements by 

, 3^]. During each deposition run, up to 8 pieces of the substrate measuring 



Rutherford backscattering C35L 
approximately 5x5 mm 2 were placed such that a number of different irradiation experiments could be 
performed in vacuum. The thickness variation between the different pieces in a single run was measured to 
be less that 10%. For every film thickness we also made atomic force measurements (AFM) of the surface 
roughness and an upper limit of 0.1 db .03 nm was established for the average rms roughness over the entire 
film thickness range. AFM was performed in ambient conditions using a Molecular Imaging PicoScan AFM 
with a tip of spring constant 0.2 N/m. The measured lateral (x — y) resolution was ~ 10 nm and z-height 
resolution was ~ 0.2 nm. A Hitachi S-4500 field-emission microscope operating at 15 kV was used for 
EDS and SEM measurements. The spatial resolution of this microscope was ~ 1.5 nm. 

The films were irradiated in vacuum at normal incidence by a laser beam of variable beam size controlled 
by a 500 mm biconvex lens, with a maximum size of 3 x 3 mm 2 . The laser used was a Spectra Physics 
Injection seeded Lab- 130-50 Nd:YAG laser operating at its 4 th harmonic with wavelength A = 266 nm at 
50 Hz repetition rate. Based on the thermal model (Sec. IIIIBI ). the metal film is expected to return to room 
temperature within approximately 500 ns and therefore each laser pulse could be viewed as being thermally 
independent of one another. The laser beam was plane polarized with a coherence length of ~ 2 m and had 
a Gaussian temporal shape. The pulse width t p , as measured by the full-width-half-maxima of the temporal 
profile, was 9 ns. The spatial profile of the beam was a circular multimode energy distribution established 
by the diffraction output coupler of this laser. The spatial energy distribution was best described by fitting 
to a Gaussian shape which gave a 85% fit to the centroid of the distribution in each direction with a full 
width at half maxima of 4 mm. In the studies reported here, all morphology investigations were reported 
from an approximately 50 x 100 am 2 region of the irradiated area within which the laser energy density 
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variation was < 10% in each of the two dimensions. 

Two types of laser irradiation experiments were performed. The first was determination of the melt 
threshold laser energy density, E m . We defined E m as the minimum energy density at which a morphology 
change could be detected after the shortest irradiation time of 10 laser pulses. We also verified that for E < 
E m no morphology change could be observed even after the longest irradiation time of 10500 pulses, thus 
ruling out contributions from surface diffusion. The morphology change observed was a typical dewetting 
pattern observed in the SEM, as shown in Sec. HUB I Such an approach has been successfully used in the past 



to detect melt threshold of thin films under laser irradiation 0711 . The laser energy density was controlled 
by varying the lens-sample position and/or the laser energy, which was measured using an Ophir electronics 
PE30 A-P high-power laser detector head with an ORION power meter. In this manner energy densities 
ranging from 10 — 1000 mJ/cm? can be achieved. The second study was of dewetting pattern evolution 
as a function of various energy densities E > E m and over a range of irradiation times, as measured by 
the number of pulses n. The irradiation time typically ranged from 10 < n < 10500. The pulse number 
was controlled by a shutter in the beam path. Each independent irradiation i.e., with a different E and/or 
n were performed on different pieces of the substrate and each was characterized using SEM and/or AFM 
measurements. While the lower limit on the choice of E was set by E m , the upper limit was set by ensuring 
that after the longest irradiation times no substantial evaporation (< 15%) had occurred. This was ensured 
by measuring the concentration of the metal after irradiation via EDS in the SEM and comparing with the 
thickness calibrated values. In this manner detailed information on pattern formation was obtained as a 
function of E and n for each h under conditions not influenced by evaporation. 

B. Thermal transport modelling 

Given the applicability of continuum-level modeling established in this work, we point out that the 
nanoscopic length scales allow for certain simplifications, which in turn allowed for obtaining insightful 
analytical solutions for the film temperature as a function of time. The simplifications include: (i) h < a" 1 
and (ii) the heat loss occurred primarily through conduction through the substrate. Typically for pulse du- 
rations > 1 ns, Fourier heating theory is found to be applicable and consequently the laser heating may 
be modeled by a position-dependent source term in the transient energy equation 113811 . The heat transport 
model was developed to simulate the experiment of heating by a pulsed laser beam incident normal to the 
film-substrate bilayer with a spatially uniform energy profile in the plane (x — y) of the bilayer that was 
infinitely wide in the x and y directions. For these conditions, the relevant heat diffusion was normal to the 
film surface in the ^-direction where z = —h corresponds to the vacuum-film interface and z = the film- 
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substrate interface. The model simulated a 1-D geometry consisting of a thin, planar and continuous metallic 
layer on a thermally insulating substrate with the schematic shown in Fig. [2a). The following assumptions 
were made in formulating the heat transfer problem: (1) The substrate layer was viewed as a semi-infinite 
medium (0 < z < oo). In the case of Si02 for UV irradiation, the substrate served strictly as a conducting 
medium since it was transparent to UV, and (2) The thermal contact resistance at the film-substrate interface 
was neglected which is justified aposteriori as follows. The numerically predicted thermal flux at the film 
substrate interface is of order ~ l0 10 W/m 2 for films irradiated with laser energy densities in this work 
ranging from ~ 50 — 350 mJ /cm 2 . Typical values of thermal resistance for thin films range from 10~ 5 to 
10~ n Km 2 /W 13911 . implying that an addition of thermal impedance would lead to model predictions for 
the melting threshold which are significantly lower that its experimentally measured values. 
In the above framework, the film temperature T*(z,t) is described by: 

dT* d 2 T* 

{pC e ff) m -QjT = km ~Q^2~ + ' a m eyLp(-a m (z + h)) (1) 

where the subscripts m and s denote the metal and substrate respectively, p is the density, C e // is the 
effective heat capacity, S* is the rate of energy flux absorbed from the laser into the film, f(t) denotes the 

temporal energy distribution of the laser source and a m is the absorptivity of the metal film at 266 nm. 

_(t-t p ) 2 

Specifically, f(t) = 1, 1 — H(t — t p ) and e 2-<t 2 for uniform, square-shaped and Gaussian pulses 
respectively with H denoting the Heaviside step function. We are interested in the average film temperature 
especially since the spatial variations in T* are expected to be small since h < a" 1 and thermal diffusion 
time scale < t p . An equation for T = r j\ T*dz can be obtained by integrating Eq. Q] with respect to z 
from z = — hto z = and dividing throughout by h: 

(pC' e// ) m — = -— H — (1 - exp(-a m h)) (2) 

where q s (t) = —k m (dT* /dz)o = —k s (dT s /dz)o is the conductive heat transfer to the substrate where T s 
denotes the substrate temperature field. Note that the radiative heat transfer from the film surface to the 
vacuum is typically much smaller than the conductive transport to the substrate. Hence, (dT* /dz) z= -h is 
set to 0. Thus, the lumped parameter equation is valid for all h given that T is interpreted as the average 
temperature in the film. Evidently, for h < a~} , the gradients would be sufficiently small so that this 

l-cxp(-q m 
(pC e ff)mh 



would be a good approximation to the uniform film temperature. Letting S = - Q ^ xp ^\ the time- 



dependent temperature in the metallic monolayer can be expressed as: 



§ = 3 . m -*#L^ (3, 
at (pC e ff) m h 
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where S is the rate of energy absorbed from the laser, averaged over the film. For a Gaussian shaped pulse 
we have: 

s *=( (1 - fl(ft)) -vf^) (4) 

and for constant heating or a square pulse 

S*=(0-R(h))-C~^ (5) 
where R is the /i-dependant reflectivity for the laser at normal incidence. The general expression for the 



effective R value for an absorbing film on a transparent substrate has been derived previously H40TI . Fig. 
[TJb) presents the change in R for Co-Si02 over the film thickness range of interest, and could be well 
represented by a saturating exponential of type r Q (l — e~ arh ), with r a = 0.44 and a" 1 ~ 15.5 nm. Further, 
in the above expressions, E a is the laser energy density, a is the standard deviation chosen so that the length 
of the Gaussian full-width half-maximum is equal to the laser pulse width t p , £ is a renormalization factor 
to ensure that the total energy absorbed per pulse is identical for the square and Gaussian models. The 
temperature field in the substrate was described by the variable Q s = T s — Tq, where To = T s (z — > oo): 

d& 1 = a ^e. 

dt s dz 2 

where the substrate thermal diffusivity a s = jjtfj^ with initial/boundary conditions 

Q(z, 0) = 8(oo, t) = and — &s^f (0, t) = q s (t), which couples the substrate and film temperature 

fields. By employing the method of Laplace's transforms, the above set of equations can be solved to yield: 

T(f) = T +~(j^j Vi+j^-(e K2t •erfc(iTVt) - l^-sJ^e K2 ^-eTfc(KVt^i)-H(u-t p ydu 

(V) 

where K = ^c^//^ an< ^ t ' ie i nte £ ra l m ^1- El ma y be solved numerically. Now for the case of the 
Gaussian model, the solution for T(t) is given by: 

T{t) = T a + S- [ exp (-/ 2 • (t - uf + g ■ (t - u) + K 2 ■ u) ■ erfc(K^H) • du (8) 

where the integral in Eq. [8] may be solved numerically using standard software packages. For both the 
Gaussian and square pulse solutions, we used the Riemann midpoint rule to evaluate the integral with 
Maple™ v9.1. We note that for a uniform pulse (f(t) = 1), a simpler analytical solution of the following 
form can be obtained: 

T(t) = To + j^-(^)V-t + ^-(e^. erfc(K^) - l) (9) 
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which for t > 1/ K 2 (where 1 / K 2 ranges from 4.8 • 10~ 3 ns to 5.9 • 10 -1 ns corresponding to Co film thick- 
ness of 1 to 11 nm) would predict a temperature rise proportional to y/t. While this solution is insightful, 
it over predicts the temperature rise within the film. Hence, a more realistic Gaussian temporal distribution 
has been employed in the subsequent analysis. 

In order to validate the assumptions made in the analytic treatment, a two-phase transient heat trans- 
fer model was also solved numerically using the finite element method by using the commercial software 
FEMLab™ v3.1. Here, the assumption of homogeneous temperature in the film was relaxed. In order 
to mimic the semi-infinite substrate phase, the boundary condition T s (z,t) \ Z =5-L th = T Q was imposed 



where L t h is the thermal diffusion length of the substrate « y/t p a s . The simulations were performed 
for different spatial and and temporal resolutions to ensure numerical convergence. For temperature in- 
dependent thermophysical parameters, we observed good agreement between the analytical and numerical 
results for T within 0.5 K (see Fig. 0a)). First order phase transition was modeled by using an effec- 



tive heat capacity technique based on previous approaches 141 



4211 . namely, C e ff = C p + where 

i Te forT m -e<T<T m + e\ 
F(T) = ^ >, L is the latent heat of transformation for the film and T m 

otherwise I 

is the melting temperature of the film. The functional form of F(T) was chosen so that C e // models a first 

order phase change i.e., when |y is integrated from T m — e to T m + e the result equals unity. Here, e is a 
small temperature parameter whose value was chosen so that no oscillation in T would be observed during 
the melting and freezing transitions. A value of e = 15 K was sufficient to ensure oscillation-free results 
in the simulations reported here. The temperature dependent and independent materials parameters are tab- 
ulated in Table J) All simulations predicted the maximum temperature difference between the film-surface 
and the film-substrate interface to be <C IK. This indicated that the z-averaged temperature field obtained 
from the analytical model can be interpreted as a nearly uniform film temperature. 



DDL RESULTS AND DISCUSSION 
A. Laser energy threshold E m to observe dewetting vs. thickness 

1. Experiment 

The morphology of the films was first investigated as a function of increasing laser energy density over 
a practical time scale ranging from a few laser pulses to a maximum of 10, 500 laser pulses. The first 
observation was that below a critical E no perceptible morphology change could be observed even after the 
longest irradiation. On the other hand, above this energy, significant morphology changes could be observed 



10 



even after the shortest irradiation time (n = 10). This sharp threshold behavior implied that a phase-change 
process, such as melting, rather than surface diffusion was responsible for the morphology change. Using 
this approach, a measurable critical E was obtained for a range of film thickness. The experimentally 
determined behavior of E m vs. h is plotted in Fig. [2{a) (filled circles). This observed trend i.e., a significant 
increase in E m with decreasing h for the ultrathin films (h < 5 nm) is consistent with the prediction of 1I37I 
and the underlying mechanisms that contribute to this behavior are discussed below. 



2. Modeling 

In order to determine whether the energy E m closely correlated with the melt threshold, we used the 
analytical and numerical models to estimate the melt threshold energy. This was obtained by simulating 
the thermal behavior under the pulse for various energies. Fig. 12b) shows a plot of Eq. [8] versus time 
for various Co film thicknesses h with a laser energy density value of E = W0mJ/cm 2 estimated using 
the analytic model. The points represent Eq. [8] solved at various times and the lines represent the cor- 
responding finite element simulations. Predictions of a numerical simulation incorporating phase change 
and temperature-dependent material parameters are shown in Fig. Elc). In addition, from Eq. [U the melt 
threshold energy E m , which is the value of E needed to bring the Co film to melting temperature was cal- 
culated. The time t m where Eq. preaches its maximum temperature was first determined for the various 
film thicknesses. Then, Eq. [8] was set equal to T m at t = t m and solved for the appropriate value of E by 
successive approximations. Analytical predictions for the melt threshold energy versus h is shown in Fig. 
|2{a) as the solid line along with the numerically predicted laser melt threshold energy density (dashed line). 
From Fig. 0a), an increasing E m with decreasing h is obtained for experiment and modeling and further- 
more, the estimated E m is in excellent agreement with experiment. Based on this result, we interpreted the 
experimentally measured E m to be the melt threshold energy and all dewetting behavior was investigated 
for energies at or above the E m for the various h. 

Eq. 9 can be used to explain the intriguing thermal behavior. Specifically, the temperature rise is 
proportional to S/K oc f\(h) /2(h) where fi{h) = 1 — exp(— a m h) and /2(h) = 1 — tq + tq exp(— a r h) 
represent the contributions to the source term (S in Eq. [3]> from finite laser absorption depth of the metal and 
metal/substrate bilayer reflectivity respectively. It is easily verified that f\ increases while /2 decreases with 
increasing h. However, since a m > a r , the net effect is an increase in S/K with increasing h. Consequently, 
thicker films will experience faster heating rates. However, note that if one disregards the latter effect (of 
reflectivity), the model would over predict the temperature rise by approximately 6 and 26 percent for h = 1 
nm and 10 nm respectively. This suggests that appropriately accounting for the thickness dependence of the 
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reflectivity, overlooked in previous work 114311 . is important to arrive at predictions which are quantitatively 
consistent with experimental measurements. Moreover, the sensitivity of S/K to laser-metal interaction 
parameters implies that the thermal behavior and hence the pattern formation itself could be highly system 
dependent. As shown in Fig. QIb), the thickness-dependent reflectivity of the metal-substrate bilayer varied 



from the value for bulk S1O2 of R ~ for h m = to the value for bulk Co of R = 0.44 144J] and in the 
thickness regime of interest in this work, i.e. h m < a~}, the calculated value is significantly smaller than 
that of bulk Co. 

Another important result obtained from the thermal modeling is the prediction of an intrinsic thermal 
gradient dT/dh that can result in temperature variations along the film whenever the irradiated film has 
height fluctuations, i.e. employing a local approximation, dT/dx = (dT/dh) dh/dx. The length scale 
of the interface height fluctuations (or the interface roug hness) ~ ^Jk B T/^ 0.2 nm for a liquid Co film 
at 1700 K. While such perturbations are randomly distributed, TFH theory predicts that a wavelength with 
fastest growing timescale (or eigenvalue, see Sec. IIII B 3l below) will be preferentially amplified and become 
the dominant, observable length scale. Temperature variations could result in surface tension gradients along 
the interface which in turn will induce thermocapillary flow that can reduce (dT/ dh > 0) or enhance (dT/ dh 
< 0) the dewetting rate via thermocapillary flows that promote fluid motion from the crest (trough) to the 

trough (crest) of the perturbed interface. Specifically, dT/dh can be estimated from Eq. |9]and the resulting 

dT 

wn in l ig. rjj wucic lug uiuniaic 

(pC e ff) m 1 ~ (pC)mtp 

the rate of change in temperature, evaluated for different time-to-melt of the films, is a strong function of 
film thickness. One can clearly note that dT(t, h)/dh becomes negative at a certain critical thickness h c 
when the film melts over reasonable time scales. h c is observed to lie between ~ 9 — 13 nm corresponding 
to melting times between 1 — 9 ns. The time scale over which dT/dh will be effective is the thermal 
diffusion time scale th in the Co metal films and this time should be compared to pulse time t p . Clearly if 
TH/tp < 1, then thermal diffusion would smooth out dT/dh effects over the time of heating and its role on 
pattern formation would be negligible. As we show in Sec. IIII B 3 1 the metal dewetting is most likely due to 
a TFH hydrodynamic instability similar to spinodal dewetting. According to this this model, the dewetting 
length length scale A varies as h 2 and thus the relevant thermal diffusion time is tjj oc A 2 /D^° oc h 4 , 
where = 2.7 x W~ 5 m 2 /s is the thermal diffusivity of Co. Using the quantitative expression for 
A obtained in Sec. IIII B 31 we determined that up to thickness of h ~ 6 nm th < t p and therefore the 
intrinsic thermocapillary effects will only have a minimal influence of pattern formation for the thickness 
regime investigated here 1 < h < 7 nm. On the other hand, since th oc h 2 significant modifications to 
pattern formation for thick films is likely, including changes to patterning length and time scales. However, 



dependence on h is shown in Fig. [3j where the ordinate ^ is scaled by the numerical factor (A 2 /P)10 9 , 
where A = V^.P^ff k }± an( j p _ t are material and laser related parameters. As seen from Fig. [3j 
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a detailed investigation for thickness above 7 nm is beyond the scope of the present work. 



B. Dewetting pattern evolution with n and E 

A detailed investigation of the dewetting morphology was performed for various laser energies as a func- 
tion of the laser irradiation time n. For films in the range of 3 < h < 7 nm the patterns typically consisted 
of discrete holes at the early stages of irradiation, followed by cellular patterns at later stages and eventu- 
ally nanoparticles which remained stable to irradiation. On the other hand, for the films with h < 3 nm, 
the morphology consisted of discrete holes followed by a bi-continuous structure with a final state again 
characterized by nanoparticles. It is widely accepted that the dewetting morphology can arise from three 



mechanisms II45I1 . (i) Homogeneous nucleation and growth in which holes appear randomly in location and 
time on the surface. Therefore, no characteristic length is present in this type of dewetting 122J]. (ii) Hetero- 
geneous nucleation and growth due to defects, impurities or other experimentally imposed heterogeneities. 
Here, a characteristic length scale could appear at the early stages of dewetting due to ordered nucleation 
sites. For instance, in ion-irradiation induced dewetting, the average molten zone of an ion imposes a char- 
acteristic length scale in dewetting J22I. (iii) TFH instabilities such as the one associated with the dewetting 
of spinodally unstable systems. The resulting patterns were characterized by a well-defined length scale in 



the hole spacing and/or size 12111 . Since we have ruled out spatially ordered heterogeneities on the surface 
and in the film microstructure, the results presented below were analyzed to distinguish between dewetting 
by homogeneous nucleation and TFH instabilities that establish a characteristic length scale. In order to do 
this, the primary characterization was evaluation of the patterning length scales obtained from the power 
spectrum (PS) evaluated by performing fast fourier transforms (FFT) of the morphology of the dewetting 
patterns. 



1. Dewetting morphology for h < 3 nm 

Fig. ffia) to (d) capture the the morphology of a 2 nm thick Co film as a function of increasing number of 
laser pulses at energy 200 mJ /cm 2 . Discrete holes were visible after the shortest time (Fig. @la)) with the 
patterns changing to a bi-continuous structure (Fig. @fc and d) and finally into nanoparticles (Fig. HJd). A 
comparison of the density of features in Fig. @Ic) an d (d) indicates that the nanoparticles form by change of 
spaghetti-like features into a single particle rather than break-up of the spaghetti-like features into multiple 
particles, such as in a Rayleigh-like process. As shown in Fig. [Ha), the nanoparticles had a monomodal 
size distribution. The spatial behavior of the patterns could be understood from the PS, as shown in Fig. 
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[Ua) to (d). During all observed stages of the pattern evolution an annular shaped PS was visible, suggesting 
that a band of characteristic spatial frequencies was present in the patterns. Qualitatively, the radius of 
the annular region changed with increasing irradiation time reflecting the change in morphology of the 
dewetting patterns. Fig. |3a) to (d) show the pattern morphology after 100 pulses as a function of various 
energies for the 2 nm film. The general characteristics of the morphology were similar to those observed as 
a function of n, as shown in Fig. 0]and the primary role of increasing energy appears to be an enhanced rate 
of achieving the nanoparticle state. 

2. Dewetting morphology for 3 < h < 7 nm 

For the thicker films the intermediate stages of dewetting differed considerably in their morphology 
from Fig. @]but still retained a characteristic spatial frequency. Fig. [6fa) to (d) capture the morphology of 
a 4.4 nm thick Co film as a function of increasing number of laser pulses at energy 93 mJ /cm 2 . Discrete 
holes were visible after the shortest time (Fig. a)) with the patterns changing to a cellular structure (Fig. 
0c)) as the number of holes increases. Continued irradiation results in the metal being pushed to the edge 
of the holes and the patterns consisted of large polygonal structures with evidence for particle formation 
occurring preferentially at the vertices. This is more evident in Fig. 0c). After long irradiation times, 
a stable nanostructure morphology was observed (Fig. 0d)). This behavior is consistent with particle 
formation via capillary draining of the liquid bridges into the vertices (Sec. MI B 3I >. 

Analysis of the PS for each pattern again showed an annular structure indicating spatial order in the 
dewetting patterns, as shown in Fig. 0a) to (d). The radial distribution function g(k) for each of these PS 
and the curve fit used to obtain L are shown in Fig. Eta). The cumulative behavior of L vs. n for this film 
is shown in Fig. Hb). As in the case of the 2 nm film, L increases as the holes appear and merge to form 
closed ring-like polygons. In this scenario, L represents the average diameter of the holes or polygons. 
However, unlike the case of 2nm film, the appearance of the nanoparticles decreases the characteristic 
length scale. This can be understood by noting that the final length scale represents the average NN spacing 
of the nanoparticles which formed preferentially at the vertices of the polygons. Since the average vertex 
spacing is smaller than the average diameter of the polygons therefore L decreased as polygons changed 
into nanoparticles. The nanoparticles once again had a monomodal size distribution, as shown in Fig. 0b). 
We also investigated the role of laser energy on the pattern formation. Similar to the case of the 2 nm film, 
we observed that the primary effect of increasing energy is an enhanced rate of achieving the nanoparticle 
state. This effect is better represented in Fig. Hb) where the characteristic L determined from the g(k) is 
plotted as a function of n for various laser energies E. As the laser energy increases, the observed value of L 
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is larger suggesting that the pattern ripening rate increased. However, R which represented the nanoparticle 
morphology after long irradiation times was independent of the laser E (Sec. IIII CI ). This result once again 
confirmed that the final nanoparticle state is robust to the laser parameters. 



3. Quantitative evidence for TFH instability 

From the above results we suggest that the pulsed laser melting initiates the nonlinear dewetting TFH 
instability and subsequent pulses evolve this instability to a well-defined final state. Further, each observed 
state is characterized by a narrow band of spatial frequencies indicating spatial order over the entire dewet- 
ting process. Since we ruled out heterogeneous nucleation from ordered surface defects, the above results 
point strongly to dewetting by an instability. A more quantitative evidence for this hypothesis can be pur- 
sued on the basis of predictions of the linear TFH theory according to which the characteristic time scales 
and length scales of dewetting are well-defined functions of h. From our previous results, the most ro- 



bust measure of the dewetting length scale appeals to be the final nanoparticle state Bill . Therefore, we 
studied the characteristics of this final nanoparticle state including, R, D, contact angle function f(9) and 
the particle areal density N as a function of film thickness h to determine pattern formation via the TFH 
instability. 

1. Rvsh 

According to the linear TFH instability, the characteristic spinodal length A as a function of film 



thickness can be expressed as 

113 



Hh) = \l^ 1 h 2 = Bh 2 (10, 



where 7 is the surface tension of the metal and A is the Hamaker constant. Since we interpret the 
final nanoparticle morphology as resulting from the dewetting instability, we assigned the observed 
NN spacing R to be proportional to A as R = aA, where a is the proportionality factor. This is 
reasonable given that the average nanoparticle spacing is related to the average size and/or spacing 
between the ordered features arising from the dewetting process. We determined the behavior of 
R as a function of film thickness, as shown Fig. Ha). The power law fit to the spacing R yielded 
a behavior expressed by R(nm) = 25.7/i 198±0 3 with the film height expressed in nm. The NN 
spacings varied from ~ 40 to 1000 nm for the films studied. The exponent of the film thickness of 
1.98±0.3 is in excellent agreement to the theoretical value of 2. From this experimentally established 
pre-exponential factor of 25.7 nm -1 the Hamaker constant was estimated to be A = 1.41 x 10~ 18 J 



which is of the right order of magnitude [46]. 
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2. Dewetting instability time scale td 

In the multiple pulse laser melting dewetting scenario, the morphology change occurs over multiple 
cycles of phase change, i.e. melting and resolidification. Clearly, once the instability is initiated, in 
order for it to be fostered by the multiple pulses, the time scale of dewetting T£> must be larger than 
the lifetime of the liquid tl in each pulse. From the linear TFH theory, m can be expressed as: 

where n = 4.46 x 1CT 3 Pa — s is the viscosity of liquid Co at its melting point. Using the experi- 
mentally determined value of A we find that T£> ~ 25 ns for a 1 nm thick film and increases as h 5 . 
From our modeling, we determined that at the melt threshold of ~ 338 mJ/ cm 2 for the 1 nm film, 
the liquid lifetime is ~ 4 ns. The typical liquid lifetime was obtained from thermal profiles of the 
type shown in Fig. Ob) and (c). The liquid lifetime for the analytic result (Fig. 0b)) was defined 
by placing a line at T m = 1768 K, the Co melting temperature, and noting the time between both 
intercepts. For the numerical model (Fig. Etc)), the liquid life time is defined as the time needed for 
the system to progress from the onset of melting to the completion of freezing. From this analysis, 
the liquid lifetime for the range of films studied here varied between 2 — 15 ns, as determined from 
the thermal model, and clearly is smaller than the dewetting time estimated from Eq. [TTJ 

We also verified from the thermal modeling that the heating and cooling rates were extremely rapid 
and since the time between pulses was large (20 ms), each pulse was thermally independent of the 
other. To extract the rates, first, plots of the type show in Fig. Ob) were generated for Co film thick- 
ness of 1, 2, 3, 5, and 11 nm with laser energy density values chosen so that the peak temperature 
would be ~ 2000 ± 100 K. The heating rate was estimated from the time At to reach maximum 
temperature T m as (T m — To)/Ai and the model predicted typical rate of temperature rise of the 
bilayer system to be ~ 150 K/ns for the metal-film thicknesses of interest. For the cooling rates, 
and a second order exponential fit to the cooling portion (e.g. Fig. [Tola)) and a typical time constant 
of ~ 130 ns was estimated for the temperature to fall by 1/e of its peak value. In addition, within 
approximately 500 ns the bilayer system reached 98% of its final value giving a time averaged cool- 
ing rate of ~4 K/ns. Clearly, each laser pulse could be treated independently of one another. These 
time scale estimates support a dewetting scenario in which the instability is initiated by a number of 
pulses and then dewetting progresses by the subsequent pulses to the final stable nanoparticle state. 

3. D and f(9) vs h 
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An analytical relationship between the diameter of a particle D and the NN spacing R can be derived 
based on volume conservation [20]. Assuming that a circular section of the film with a volume 
approximated as a cylinder of diameter R and thickness h contributes to one spherical nanoparticle 
of diameter D, then volume conservation suggests that: 

ltfh=%.f(e)£ (12) 

4 3 J w 8 

where f(9) is the geometric factor based on the particle contact angle 9 that determines the fraction 
of the sphere that lies above the substrate ( for the case of 9 = 180°, f(9)=l). Using the characteristic 
length from Eq. [l0]in Eq. [12] we get: 

D = (xSr) 1/3/l5/3 = Ch5/3 W 

In Fig. |2b) the average particle diameter D is plotted as a function of h. D varied from ~ 30 
to 250 nm for the thickness range studied here. The power-law fit to our experimental data gives 
D(nm) = 14.8/i 16±0 3 with the film height expressed in nm. The exponent of 1.6 ± 0.3 is in 
excellent agreement with the theoretical estimate of 5/3. By using the experimentally determined 
pre-exponent of 14.1 ran" 2 ' 3 in Eq. [13J we determined the geometric factor to be f(9) = 0.361 
giving an effective contact angle for the nanoparticles of 74°. Given the simplicity of the model 
relating the instability length scale to the particle diameter and the uncertainties in the values of R and 
D of approximately 30%, this experimentally extracted value of the contact angle is in reasonable 
agreement to the previously reported value of 100 ± 26° measured by atomic force and scanning 



electron microscopy 13211 . 
Areal density N vs h 

The areal density N of particles can be defined as: 



N = 1/(^A 2 ) (14) 

which using Eq. [lO]results in a thickness dependence given by: 

N(h) = ^-h~* (15) 

Plotted in Fig. |9]c) is the variation in the observed particle density with h. A power law fit gives = 
1.5 x io~ 3 /i~ 3 ' 8±a8 #/nm 2 with the film height expressed in nm. The exponent of —3.8 ± 0.8 is in 
fair agreement with the expected value of 4 (with the uncertainty primarily arising from propagating 
errors in the various quantities) while the pre-exponent is consistent with the value obtained by 
substituting the experimentally determined A and 7 in Eq. [15] 
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5. Effect of large capillary velocity on pattern formation 

From the quantitative results above, the classical TFH dewetting instability appears to be a suitable 
explanation of the pattern formation in these nanoscopic Co films in the regime of weak thermocap- 
illary effects. At this stage it is also important to emphasize the difference between the dewetting 



patterning behavior observed in the metal films reported here to those of typical polymer films M20I1 . 
Previous studies of the break-up of cellular structures made from typical polymeric liquids show that 
the sides of the polygons break into a chain of droplets by the Rayleigh instability and as a con- 
sequence the particle diameter will follow the trend D ~ /i 3 / 2 12011 . However, for the metal films 
investigated here, the final break-up of the cellular regions most likely occurs via capillary flow of the 
liquid into the vertices of the polygon. We arrived at this conclusion based on the close correlation in 
the spatial density of the vertices and particles (Fig. Oc) and (d)) or the spaghetti-like features and 



48] 



particles in Fig. @Jc) and (d). We attribute this to the extremely large capillary velocities 1471 
of the liquid metals given by the quantity v c = 7/77, which for Co is ~ 422 m/s. This implies that 
capillary draining from the liquid bridges connecting the vertices (reservoirs) which are separated 
by a distance of 1 \im (which is smaller than the largest NN spacing reported in this work) is ~ 2.5 
ns. This is smaller than the (Rayleigh) ripening time (TR ay ~ 27r(pr 3 / C7) 1 / 2 ), where p is the liquid 
density, r is the cylinder radius and c is a factor of 0(0.1 ) determined by contact angle and boundary 
conditions, for surface perturbations of the liquid bridges, which typically ranged between several 
to > 10 ns for the bridge sizes observed here 149(1 . In this respect, the final stages of breakup is 
different from that typically observed in common polymer films investigated, such as polystyrene, 
in which the cellular rims break-up via a Rayleigh-like instability Moreover, we point out that 
the Ohnsorge number, Oh = rr^p tlrat represents the ratio of viscous to capillary forces is typ- 
ically << 1 for typical O(10 nm) values of the liquid bridge diameter, D. This further underscores 
the dominance of strong capillary effects in the breakup of polygonal structures. For more viscous 
liquids (e.g. polymers) with typically smaller interfacial tension as compared to metals, the relative 
importance of viscous and capillary forces can change by orders of magnitude, leading to different 
mechanisms of breakup, e.g. via Rayleigh-like process. Another factor that can influence the particle 
size distribution is coarsening due to surface diffusion. However, as seen from Fig. |7{b) changing 
the laser energy, which changes liquid temperature and lifetime (Sec. MI CI ), did not change the final 
spacing. This was attributed to the fast thermal processing that minimizes particle coarsening. 
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C. Role of Co liquid temperature and lifetime on dewetting 

From the experimental results of dewetting presented in the previous section, the length scales L during 
the early stages of dewetting increased with increasing laser energy following similar irradiation times 
(Fig. [7]). From our thermal modeling, increasing laser energy increased the Co liquid temperature and the 
liquid lifetime. Further, the length scale of dewetting, A and the dewetting timescale td have temperature- 
dependent materials parameters (i.e. 7 and 77). Therefore the observations of Fig. |7]can arise from three 
factors: (i) from a change in the dewetting length scale A, (ii) from an increase in liquid lifetime tl per 
pulse, and (iii) from a change in the time scale for ripening tr of the dewetting structures. However, Fig. 
|7] also indicates that changing the laser energy had no influence on the final nanoparticle characteristic, 
implying that no permanent change to A must occur with change in liquid T^. In this section we show that 
the change in A for Co metal is experimentally indistinguishable over significant T changes (of up to 500 
K above T m ). From theoretical arguments we found that the increase in L with increasing E following 
irradiation by a small number of pulses n (i.e. early stages) arises primarily due to variations in tl and 
its linear behavior with E. This result was also consistent with experimental observations based on the 
measurement of L vs n for various energies and thickness. 



1. Role of temperature on A 



The T-dependence of the dewetting length A, which comes from the surface tension 7(T) and the 
Hamaker constant A(T), can be expressed as: 



A(h,T) 



'lforMT)^ 



A(T) 



(16) 



The variation in temperature can be expressed as: 



dA 
dT 



VWtt 3 



1_ _ 1/2 fry 7 1/2 dA 



2A 1 / 2 
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dT 2A 3 /2 dT 



(17) 



Eq. [17] can be simplified for metals on Si02 system by noting that the Hamaker constant for metals is 
virtually independent of temperature 1 3l and so M ~ resulting in: 



dT 



1 



r7 



1/2 #7 



h 1 



2A 1 / 2 1 dT 

For a temperature rise of AT over the melting point, the change in length can be expressed as: 

AA(AT) = Pj^-^IaT = -2.66 x 10~ 4 A(T m )Ar 
27(T m ) dT 



(18) 



(19) 
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where A(T m ) is the length at the melt temperature and "y(T m ) and d'y/dT for Co are taken from TableU 
From this we see that A decreases as Tl increases, but this decrease is small and for temperature rises of up 
to AT = 500 K is only ~ —13%. This change is smaller than the spread in the experimental values of R of 
±20% as measured from the g{k). Therefore, for conditions under which the Co liquid T rise is ~ 500 K 
above T m the change in length scale will be experimentally indistinguishable. 



2. Role of T on ripening rate L jn 



According to the linear 
to the dewetting time td Il5ll 



theory TFH theory, the pattern ripening time tr is expected to be proportional 
and so we can express the T-dependent ripening time as: 



TR(KT)= „ TD= ^HT)fT) hS (2fl) 
From this we can also define the ripening rate as a oc — , which can be evaluated using Eq. fT6land|20l 



as: 



a(h, T) = = M 7 (T)-V^ ( r)-^V2 /t -3 (21) 



where M = ^^^r is a proportionality constant with unknown numerical factors of k and u>. Using 
these expressions, the experimentally observed dewetting length L for a film of thickness h at a given laser 
energy E > E m following irradiation by n pulses can be expressed as: 

L(h, E, n) = a(h, E)t l (E, h)n (22) 

where tl, the lifetime of the liquid, is a function of the laser energy and film thickness. In order to evaluate 
the expected behavior, each of the temperature (and laser E) dependent terms on the rhs can be evaluated: 

1. Behavior of ripening rate a(h,T) for Co: 

The temperature dependence of a can be expressed analytically as: 



da(T, h) _ B ft _ 3 



Itt Vt 
2 7 T] 



(23) 



dT ^ 1 / 2 r] 

where, <y T = ^ anc j ^ T = *L an( j the temperature dependence of the Hamaker constant has been 
neglected. The behavior of <r(T) with temperature can be evaluated for Co metal by using the 
materials parameters from Table J] Fig. [TTTa) shows that Eq. [23] is well described by a linear 
behavior which increases with increasing temperature of the form: 

a(T, h) = fJA^h- 3 [0.29T - 350] = a T (h)T - a (h) (24) 
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2. Dependence on liquid T and lifetime tl on E: Thermal modeling 

From thermal profiles of the type shown in Fig. E£a), the peak temperature and liquid lifetime were 
calculated as functions of h and E. Fig. |2jb) shows thermal profiles of the film temperature for 
various film thickness for a laser energy density value of E = 125 mJ/cm 2 . The liquid lifetime for 
the analytic model was defined by placing a line at T m = 1768 K, the Co melting temperature, and 
noting the time in between both intercepts. For the numerical model, the liquid lifetime is defined as 
the time needed for the system to progress from the onset of melting to the completion of freezing. 
From Fig. [2jb) we see that introducing phase change into the model causes an appreciable change 
from the predicted temporal profile by Eq. [8] The portions of the figure where the temperature 
plateaus represent the film melting and freezing. Fig. fTOl b) presents the liquid lifetime and peak 
temperature versus laser energy density trends using the numerical and analytic models. From the 
figure, the numerical (open symbols) predictions for liquid lifetimes are about 20% lower then those 
predicted by the analytic model (closed symbols) for films greater than 1 nm. Further from Fig. 
[TOlc). the peak temperatures predicted by the numerical simulation are lower by ~ 200 K than those 
predicted by the analytic model. Both these differences arise from the fact that the phase change 
in the numerical model consumes more energy and takes more time. Fig. [TOlb) and (c) present the 
liquid lifetime and peak temperature as a function of laser energy density above the melt threshold 
for various film thicknesses. From thermal cycles of the type shown in Fig. [2j a linear relationship 
was observed between the liquid life-time and peak temperature for small changes in laser energy 
density above E m and can be expressed as: 



where the slopes A and C and intercepts B and D are thickness dependent positive quantities. Their typical 
values ranged between: 0.07 < A < 0.3 ns - cm 2 /mJ; -22 < B < -17 ns and 2.3 < C < U.5K - 
cm 2 /mJ; and D ~ 301 K for thickness between 1 < h < 11 nm with < 5% difference between the 
numerical and analytical estimates. In this laser energy density regime, while the temperature of the Co film 
is above the melting temperature, there is minimal evaporation over the relevant time scales. 

Using Eq. [25] in Eq. [22] and Eq. [26] in Eq. [24] we can express the rate of change of dewetting length L 
with n as: 



T L (E,h) = A(h)E — B{h) 



(25) 



T L (E,h)=C{h)E + D{h) 



(26) 



L(E,h) 



pE 2 +qE + r 



(27) 



n 
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where p = otC A; q = <jt(DA — BC) — a a A; and r = —B(cjtD — a a ) are the resulting temperature 
and thickness coefficients. Equation [27]can be further simplified by analysing the magnitudes of the various 
coefficients using the range of values for the coefficients A to D which showed that for films in the range 
of 1 to 1 1 nm at the melting point of Co, the coefficient p < < q and the quadratic term in Eq. [27] can be 
neglected to give: 

!^> =qE + r (28) 
n 

Therefore the experimentally observed trend in the ripening rate with laser energy should be linear. As 
shown in Fig. PTTT b). the L/n for for various film thickness can be described by linear fits (dashed lines). The 
physical significance of this linear behavior is the weak temperature dependence of the material parameters 
of Co. Therefore the primary contribution to the ripening rate due to increasing energy is the increasing 
lifetime of the liquid, which as shown in Fig. [TOl b). varies linearly with E. This result also implied that 
understanding the thermal behavior could result in predictive models describing the quenched-in dewetting 
morphologies and length-scales when such nanoscopic films are subjected to multiple instances of phase 
change. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the dewetting behavior of nanoscopic Co metal films with thickness h < a" 1 
on Si02 substrates following multiple ns pulse laser irradiation using 266 nm UV light through experi- 
ments, first-principles thermal transport modeling and analysis based on TFH theory. We determined that 
by accounting appropriately for the changes with respect to h in the reflectivity of the substrate-metal bi- 
layer, quantitative agreement with experiments can be obtained for the /j-dependence of the melting energy 
threshold predicted by continuum-level analysis. This observation corroborates atomistic simulations of 
heat transport for simpler systems under shear flow where the results were found to agree with the predic- 
tions of continuum-level analysis for liquid layers with several atomic layers. For the metal films used in 
this work, this translates to h < 1 nm. A novel effect, the development of intrinsic thermal gradients dT/dh 
whose sign is positive for h < a critical value h c rj 9 — 13 nm and negative for h > h c , was predicted by the 
thermal model. This effect arises due to the strong film thickness- and time-dependence of temperature for 
liquid metal nanolayers and fast thermal cycles. This intrinsic effect can lead to thermocapillary effects that 
could modify pattern formation. However, for the experimental thickness regime investigated, the thermal 
diffusion time scales th < pulse time t p = 9 ns and therefore the influence of dT/dh is minimal. However 
since th oc h 4 , thermocapillary mechanisms could be important for thicker films, a hypothesis that needs 
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to be verified by future experimental studies. Because of the extremely fast thermal cycles (heating/cooling 
rates of O(10 10 K/s)) and short liquid lifetimes tl of 2 — 15 ns, multiple pulses were required to observe 
the thin film instability with characteristics dewetting times td > 25 ns. The dewetting patterns showed 
short range order and typically evolved with increasing number of laser pulses from discrete holes to bi- 
continuous or cellular structures and eventually into nanoparticles. The large capillary velocity and short 
patterning length scales in the metal layer resulted in nanoparticle formation by capillary flow rather than 
via a Rayleigh-like breakup typically observed in polymer liquids. This behavior can be attributed to the 
unique thermophysical properties of the metals, including the large 7, small r\ and large Hamaker constant 
A that result in large capillary velocities and small patterning length scales. This final nanoparticle state was 
extremely robust with nearest neighbour spacing R ~ h? and a monomodal size distribution with diameter 
D ~ /1 5 / 3 respectively. This robustness was due to the relatively weak temperature dependence of the 
thermophysical parameters of Co and small particle ripening rates due to the fast thermal processing. The 
experimental ripening rate L/n increased linearly with laser energy density E at the early stages of pattern- 
ing and due to the linear increase in tl liquid lifetime with E predicted by the thermal modeling implying 
that accurate predictions and control of the quenched-in dewetting morphologies and length scales can be 
achieved in these nanoscopic systems. These results show that nanoscopic metal film dewetting following 
fast thermal processing result in novel quenched-in dewetting morphologies and length scales which can be 
accurately predicted by first-principles modeling. 
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Table Captions 

Table U Material and scaling parameters (T-dependent and T-independent) used in modeling. 
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Figure Captions 



Figure [Q (a) Schematic diagram for modeling of laser heating of a metallic nanolayer on a semi- 
infinite insulating substrate, (b) Plot of thickness-dependence of reflectivity R of Co/Si02 bilayer 
estimated from Ref. ken. 




Figure |2j Thermal behavior of the Co film on Si02 irradiated by a pulsed laser, (a) Time-dependent 
temperature obtained for the analytic and numerical model using temperature independent parame- 
ters. The symbols are from the analytical model while the lines are from the numerical calculations, 
(b) Numerical results incorporating the phase change and T-dependent parameters for various film 
thickness, (c) Laser energy density threshold for melting E m of Co films on Si02- Comparison of 
experimental measurement (solid circles) with calculations from our analytical model (solid line) 
and numerical model (dashed line). 

Figure[3j The thermal models predict an intrinsic thermal gradient dT/dh whose magnitude and sign 
was dependent on film thickness and time to melt during the ultrafast heating of nanoscopic metal 
films. As shown in the figure a critical range of thickness h c is predicted over which the dT/dh 
changes sign, and this value is dependent on the time to melt. For the thickness regime investigated 
in this work (h < 7 nm), the dT/dh dissipates in t < t p and this intrinsic effect is negligible. 

Figure |4j Dewetting pattern evolution in a 2nm Co film irradiated at 200 m J /cm 2 as a function of 
number of pulses n. The top row are SEM images while the bottom row show corresponding power 
spectrum, (a) n = 10; (b) n = 100; (c) n = 1000; (d) n = 10500. All the power spectra have an 
annular structure indicating a band of spatial frequencies implying short range spatial order. 

Figure [5J Dewetting pattern evolution in a 2 nm Co film irradiated as a function of laser energy 
density for irradiation by 100 pulses, (a) 190 mJ/cm 2 ; (b) 200 m J /cm 2 ; (c) 220 mJ/cm 2 ; (d) 
250 mJ /cm . 

Figure [6j Dewetting pattern evolution in a 4.4 nm Co film irradiated at 93 mJ/cm 2 as a function of 
number of pulses n. The top row are SEM images while the bottom row show corresponding power 
spectrum, (a) n = 10; (b) n = 100; (c) n = 1000; (d) n = 10500. Similar to the results for the 
2 nm film, all the power spectra have an annular structure indicating a band of spatial frequencies 
implying short range spatial order. 

Figure |7J Quantitative analysis of the spatial frequencies resulting from the short range order in the 
dewetting patterns for the 4.4 nm Co film, (a) The radial distribution function g(k) for the patterns 
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of Fig. [6] showing an initial increase as the holes expand and then a decrease as the nanoparticles 
start to form. The vertical lines indicate peak positions in the g[k). (b) Behavior of the characteristic 
length extracted from the peak in the g{k) as a function of irradiation time n at various laser energies. 
The spacing of the final state is independent of the laser energy. 

Figure [8j Monomodal size distribution and average diameter of nanoparticles in the final robust state 
of laser-induced dewetting. (a) An average diameter of 38 nm results for the 2 nm Co film; (b) 
Average diameter of 207 nm for the 4.4 nm Co film. 

Figure |9j Analysis of the various length scales in the final robust nanoparticle state as a function 
of film thickness, (a) The average nearest-neighbour spacing R of nanoparticles; (b) The average 
nanoparticle diameter D; (c) The average areal density of nanoparticles N in #/cm 2 . The fits show 
exponents in good agreement with linear TFH theory. The R 2 values in each graph indicate the 
correlation coefficient of the least-squares fitting. 

Figure [TOJ (a) Analytical result of the cooling portion of thermal cycle for 3nm Co on Si02 showing 
a fit using a second order exponential function. Cumulative results of the analytical (closed symbols) 
and numerical model (open symbols) using T-dependent materials parameters: (b) liquid lifetime 
vs. laser energy density; and (c) peak temperature vs. laser energy density. Closed symbols are 
analytical results while the open symbols are from numerical calculations. 

Figure [TT] Role of temperature on the ripening rate of dewetting. (a) A linear dependence of the 
spinodal ripening rate is predicted from the temperature dependence of the Co liquid material pa- 
rameters (Table. Q]). (b) The experimentally observed ripening rate L/n for various film thickness 
and energies. The symbols are experimental data while the lines are linear fits. 
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Description 


Value/Expression 
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